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Abstract-This paper presents a method for designing of 
proportional-integral-derivative (PID) controller based on 
unified power flow controller (UPFC) using chaotic harmony 
search algorithm (CHS) for damping of low frequency 
oscillations (LFO) in a power system. The main aim of this 
paper, optimization of selected gains with the time domain- 
based objective function which is solved by chaotic harmony 
search algorithm. The performance of the proposed single- 
machine infinite-bus system (SMIB) equipped with UPFC 
controller- based PID controller is evaluated under various 
disturbances and operating conditions. The effectiveness of 
the proposed UPFC controller based on PID controller to damp 
out of oscillations over a wide range of operating conditions 
and variation of system parameters is showed in simulation 
results and eigenvalue analysis. 

Index terms- UPFC, Chaotic harmony search, PID controller, 
SMIB system, LFO 

I. Introduction 

UPFC is a kind of FACTs devices that it can be used for 
voltage regulation compensation, phase shifting regulation, 
impedance compensation and reactive compensation [1]. 
UPFC has some great feature but first capability of it is 
independent controlling of real and reactive power flows in 
the transmission system [2]. Now days UPFC's utilize move 
forward to the secure and economic operation in new power 
systems. Several controlling methods for FACTs devices 
have been introduced. An assessment of the effects of power 
system stabilizer (PSS) and static phase shifter (SPS) control 
using simulated annealing (SA) in [3] has been carried out. 
Problem of the STATCOM state feedback design based on 
zero set concept coming in [4] . Particle swarm optimization 
(PSO) for achieved output feedback parameters of UPFC in 
[5] is used. 

In this article a single machine infinite bus (SMIB) as a 
power system is considered. After making a linear system 
around the operating condition with a disturbance in different 
loading situation that are converted to optimizing problem. 
For solving these kinds of problem, several algorithms are 
recommended in the above but in this paper chaotic harmony 
search (IHS) technique is used. By considering available 
parameter Au as a input of PID controller this procedure is 
carried out. In this paper two UPFC inputs applied 
independently that connected to the output c itroller. 
Finally the effectiveness of this work by result evaluation 
and comparison of performance indices is showed. 
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II. Chaotic harmony search algorithm 

This section describes the proposed chaotic harmony 
search (CHS) algorithm. First, a brief overview of the IHS is 
provided, and finally the modification procedures of the 
proposed CHS algorithm are stated. 

A. Improved Harmony Search Algorithm 

Fundamental of this method is based on the concept in 
search of a suitable state in music. In that mean's how a 
musician with different search modes to reach its desired 
state. To implement the above concepts in form of algorithm 
ahead will several steps. This process generally in the form 
of the following five steps will be implemented [7- 9]. 

1. Initialize the problem and algorithm parameters. 

2. Initialize the harmony memory. 

3. Improvise a new harmony. 

4. Update the harmony memory. 

5. Check the stopping criterion. 

A.l. Initialize the problem and algorithm parameters 

At this stage initialize all the algorithm parameters are 
carried out. These parameters can be defined by the following 
terms [7- 9]: 

HMS: Harmony memory size or the number of solution 

vectors in the harmony memory 

HMCR: Harmony memory considering rate 

PAR : Pitch adjusting rate 

NI : Number of improvisations 

N : Number of decision variables 

A.2. Initialize the Harmony Memory 

The harmony memory (HM) is a matrix include of 
independent variables which in this stage are selected 
randomly. Each row of this matrix indicates a solution to the 
problem. 



HM = 



r HMS 



CO 
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A.3. Improvise a New Harmony 

Producing a new harmony is called improvisation. The 
new harmony vector x =( x ', x 2 ',. . ., x ') is achieved by three 
parameters: harmony memory consideration rate, pitch 
adjustment rate and random selection. In the improved 
harmony search algorithm two parameters in each iteration 
are changed. These parameters are pitch adjustment rate 
(PAR) and bandwidth (bw). The form of the changing these 
parameters is shown at below equations: 
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PAR=PAR mLn + 



Where: 
PAR 
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pitch adjusting rate for each generation 
minimum pitch adjusting rate 
maximum pitch adjusting rate 
number of solution vector generations 
generation number 



bw(gn) = bw mflX . e 



(3) 



Here 

bw(gn) bandwidth for each generation 
bw minimum bandwidth 

nun 

bw maximum bandwidth 

max 

A.4. Update the Harmony Memory 

In the updating harmony memory stage, with placing new 
value ( x ' = ( x , x_ , . . . , x N ') ) on the objective function if objective 
function value is better than the previous case it must be old 
harmony out of memory and its new value will replaced. 

A. 5. Check Stopping Criterion 

If the stopping criterion (maximum number of 
improvisations) is satisfied, computation is terminated. 
Otherwise, steps 3 and 4 are repeated. 

B. Proposed Method 

In numerical analysis, sampling, decision making and 
especially heuristic optimization needs random sequences 
with some features. These features consist a long period and 
good uniformity. The nature of chaos is apparently random 
and unpredictable. Mathematically, chaos is randomness of 
a simple deterministic dynamical system and chaotic system 
may be considered as sources of randomness [7]. There are 
many methods for generating of random numbers but here 
used sinusoidal iterator for chaotic map [7]. It is represented 
by: 

X n+ =ax n 2 sin(^x n ) (4) 

When a = 2.3 and x =0.7 it has the simplified form represented 
by: 

X ,= sin (tix) (5) 

n+l v n' v ' 

It generates chaotic sequence in (0, 1). Initial HM is generated 
by iterating the selected chaotic maps until reaching to the 
HMS as shown in flowchart Fig.l [7]: 
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Figure 1. Flowchart of the chaotic harmony search 

III. Modeling the Power System Under Study 

The system in this article is a single machine infinite bus 
(SMIB) with UPFC installed as shown in Fig. 2 [5]. For 
analysis and enhancing small signal stability of power system 
through UPFC, dynamic relation of system is needed. To have 
these relations through park transformation and ignoring 
resistance of both boosting and exciting transformers can be 
achieve following equations [5, 6]: 



/■ V Etd1 



tqj \X E JWJ 



+ 



+ 



niEcasCaE^'dc 



mBsln(fiB]vd r 



(6) 



0) 



4C dc 



3m B 



4C 



SLtlS 
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where v , i , v , and L are the excitation voltage, excitation 

EC E' Be B o ' 

current, boosting voltage, and boosting current, respectively; 
C dc and v Jc are the dc link capacitance and voltage [5]. The 
non linear model of the SMIB system introduced Fig. 2 as 
shown following equations: 



lb — (i) D (di — 1J 



17 



(9) 
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Figure 2. SMIB power system equipped with UPFC 

A. Power System Linearized Model 

With applying linearization process around operating 
point on under study system, state space model of system 
according below relation can be achieved. 

£ = A3£+BU (13) 
Where, the state vector x, input vector u, A and B are: 
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IV. PID Controller for UPFC 

The eigenvalues of the state matrix A that are called the 
system modes define the stability of the system. By using 
PID controller, can move the unstable mode to the left-hand 
side of the complex plane in the area of the negative real 
parts. A PID controller has the following structure: 



Gpm — Kt, + ^ + K d s 



(IS) 



By properly choosing the PID gains (K , K ., K d ), the 
eigenvalues of system are moved to the left-hand side of the 
complex plane and the desired performance of controller can 
be achieved. UPFC controllers (m B , m E , 5 S and 5 £ ) are four 
parameters that can help them, reach the above goals 
mentioned but in this paper the two UPFC s controller inputs 
(5 E , m B ) separately under various operating conditions is 
considered. In this study for solving optimization problem 
and reach the global optimal value of coefficients K, CHS 
algorithm is used. 



Figure 3. UPFC with PID controller 

in form Integral of Time Multiplied Absolute value of the 
Error (ITAE) as the objective function for the algorithm is 
selected. The objective function is defined as follows: 



Objective function = / D t|Adj|dt 



(19) 



In the above equations, t . is the time range of simulation. 
The design problem can be formulated as the following 
constrained optimization problem, where the constraints are 
the controller parameters bounds: 



<K < K% 



fif* <Ki< 
Kg* <K d < Kg™ 



(21) 
(22) 
(23) 



Typical ranges of the optimized parameters are [0.01-150] for 
K p and [-5-5] for K. and [-10- 10] for K d . The optimization of 
UPFC controller parameters is carried out by evaluating the 
objective cost function as given in Eq. (19). The operating 
conditions are considered as, Normal, light and heavy load 
respectively: 

P = 0.80 pu , Q = 0. 1 14 pu and x L = 0.5 pu. 

P = 0.2 pu, Q = 0.007 pu and x L = 0.5 pu. 

P= 1.20 pu, Q = 0.264 pu and x L = 0.5 pu. 

In order to acquire better performance, the parameters of CHS 

algorithm are showed in table i. 

Table i. 

Value of CHS Algorithm Constant Parameters 





HM5 


HMCR 


PAR„ 


PAR_ 






valu ? 


20 


0.9 


0.25 


0.55 


le-4 





In optimizing values of the controller parameters, algorithm 
must be repeated several times. Finally values for PID control 
gains are selected. Optimal values in normal load are shown 
in the table II: 
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Table ii. 

The Optimal Gain Adjusting of The Proposed Controllers 



controller 




Ki 






95.244 


0.B45 


1.024 




145.81 


1.027 


-9.680 



V. Time Domain Simulation 

System performance with the values obtained for the 
optimal PID by applying a disturbance at t = 1 for 6 cycles is 
evaluated. The speed and generator power deviation at normal, 
light and heavy load with the proposed controller based on 
the 5 E and m B are shown in Fig. 4 and Fig. 5 respectively. To 
demonstrate performance robustness of the proposed 
method, two performance indices: the ITAE and Figure of 
Demerit (FD) based on the system performance characteristics 
are defined as [8]: 

ITAE = 10,000 / D LD t(|Aw| + 0,03|APj)dt (24) 

FD = (1000 xOS] : + (3000 x US] : + T 3 : (25) 

Where (Aco) is the speed deviation, (aP ) is the generator 
power deviation, Overshoot (OS), Undershoot (US) and T is 
the settling time of speed deviation for the machine is 
considered for evaluation of the ITAE and FD indices. It is 
worth mentioning that the lower the value of these indices is, 
the better the system response in terms of time-domain 
characteristics. It can be seen that the values of these system 
performance characteristics with the PID § E based controller 
are much smaller compared to PID m B based controller. 
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Figure 4. Dynamic responses for speed deviation with normal (a), 
heavy (b) and light load (c): solid (8 E based controller), dotted (m B 
based controller), dashed (without controller) 
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Table hi. 

Eigenvalues of System in Different Operating Conditions for 3 e based 
controller 



Normal 
P = 0.B (pu) 


Light 
P = 0.2 (pu) 


Haavy 
P=12 (pu) 


-19.2472 


-19 212 


-19.2737 


-2.205 = 3. B9j 


-1.BB6 =4.252j 


-2.56E = 3.226j 


-1 2629 


-1.277 


-1.3248 


-0.064 


-0.0429 


-0.0584 



Table iv. 

Eigenvalues of System in Different Operating Conditions for m E based 

CONTROLLER 



Normal 
P = 0.B (pu) 


Light 
P = 0.2 (pu) 


Haaw 
P=1.2 (pu) 


-19 323 


-19.233 


-1.9.377 


45.629 = 4.372] 


45.459 =4.580] 


45.748 =4. 00 Bj 
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Figure 5. Dynamic responses for generator power deviation with 
normal (a), heavy (b) and light load (c): solid (8 E based controller), 
dotted ( m B based controller), dashed (without controller) 

Table v. 

Values of Performance Indices ITAE and FD 



Control! ar 


normal 


light 


haavy 


ED 


ITAE 


FD 


ITAE 


ED 


ITAE 


& 


7.85 


8.869 


9.895 


10.07 


5.463 


7.825 




41.4 


26.06 


69.45 


3191 


29.21 


21.07 
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Conclusions 
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The chaotic harmony search algorithm was successfully 
used for the modeling of UPFC based PID damping controller. 
For each of the control signals from available state variable 
Acq is used. The efficient of the proposed UPFC controller for 
improving dynamic stability performance of a power system 
are illustrated by applying disturbances under different 
operating points. Results from time domain simulation shows 
that the oscillations of synchronous machines can be easily 
damped for power systems by this method. To analyze 
performance of UPFC's controller two indicators were used. 
These indices in terms of TTAE' and 'FD' are introduced that 
this indices demonstrate that PID with 5 E based damping 
controller is superior to m B based damping controller. 

Appendix A 

Table.vi. System Parameters 



G*n=rator 


M = B MJ/MVA, Td^ 5.044s 
Xd = Ipu. Xr = 0.6pu . X > d = 0.3pu 


Excitation system 


Ea=io, Ta=o.o; 5 


Trans formars 


X=X I =X 1 = 0.lpn 


Transmission lin? 




Op aratinE condition 


P = 0.Spu : V = =l pu 


DC link parameter 


V7>:=1 pu. Civ;=l pu 


UPFCparam-tjr 


fe=-85.3i : im = -:S2\ .T, = 0.0i 
nn = 0.4,mB. = 0.0S ; K, = 1 
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